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We investigate the existence of spiral ordering in the planar 
spin orientation of skyrmions localised on a face centered rect- 
angular lattice (FCRL). We use the non-linear sigma model 
(NLSM) to numerically calculate the minimum energy config- 
urations of this lattice around the v = 1 quantum Hall ground 
state. Our variational ansatz contains an angle 9, characteris- 
ing the FCRL and an angle q, characterising the orientational 
order. As v is increased towards one, there is a smooth transi- 
tion from the triangular lattice (TL) characterised by (9,q) — 
(120°, 120°) to FCRLs with spiral orientational order. The 
novel feature we find is that these phases are characterised by 
(9,q) values such that 8 + q — 240° (same as the TL phase). 
As v increases further towards one, there is a sharp transition 
from the FCRLs to the square lattice (SL) characterised by 
(9, q)— (90°, 180°). Consequently, the parameter 6 + q jumps 
sharply at the FCRL-SL transition and can serve as an order 
parameter to characterise it. 

The lowest energy charged excitations about the v = 1 
ferromagnetic quantum Hall ground state are skyrmions. 
They are topological objects in which the spin grad- 
ually twists over an extended region. Their spin is 
greater than 1/2 and they carry an electric charge of 
±e. Skyrmions have a topological charge which equals 
their electric charge at v = 1 [[[]-[§. Skyrmionic excita- 
tions are favoured over single particle excitations when 
the Lande g- factor is small i.e g — > Jl],||||. Skyrmions 
which are produced around the v = 1 ground state 
have been experimentally seen by in OPNMR and op- 
tical magneto-absorption experiments on Ga nuclei in an 
electron doped multiple quantum well structure JlT| , ^2| . 
They see a sharp fall in the spin polarisation of the 2D 
electron gas on either side of v = 1 , which indicates that 
a number of spins are being flipped by the addition or 
removal of a single electron, in contrast to single particle 
excitations where no extra spin is flipped. Low g exper- 
iments observe large skyrmions with a large number of 
flipped spins p3lE4[. 

The low temperature ground state of a system of in- 
teracting skyrmions is expected to be a crystalline lattice 
|^-^|. Earlier calculations performed in the mean field 
limit by Brey et. al. || suggest that the ground state 
of a dilute system of a skyrmions is a square Skyrme 
crystal with an anti-ferromagnetic order in the planar 
component of the spin of the skyrmion. Later studies 
by Green et. al. , Rao et. al. and Timm et. al. 
analyse the classical ground state of the skyrmion lat- 
tice as v — > 1 using the non-linear sigma model (NLSM) 
. Green et. al conclude that the ground state is 
a triangular lattice with Neel ordering [Q. Rao et. al. 
report sharp triangle — > square — > triangle transitions 



in the skyrmion lattice as a function of the filling fac- 
tor near v = 1 and at T=0 ||. Timm et. al. study a 
system of well separated skyrmions as described by an 
anti-ferromagnetic XY model and propose a T=0 phase 
diagram in which the TL and SL phases are separated 
by Neel ordered, centered rectangular phases ||. The 
classical and quantum phase transitions occurring in the 
Skyrme crystal has been studied using Hartree-Fock cal- 
culations by Cote et. al. fiofl . 

In this paper, we consider spiral ground states as candi- 
date minimum energy solutions of a system of skyrmions 
described by an effective classical 0(3) NLSM. Spiral 
ground states are known to be the ground states of frus- 
trated anti-ferromagnets. For example, the triangular 
lattice with ABC sublattice spin ordering is a frustrated 
spin system in which the spin at every vertex is rotated 
through 120° relative to the others. 

Our calculations are is valid for a system of overlap- 
ping skyrmions at small g. The Coulomb and Zeeman 
terms in the NLSM compete to generate a size for the 
skyrmion and to drive the TL-SL transition through the 
FCRLs. This is in contrast to the regime of Timm et. 
al. where the Zeeman contribution is neglected in the 
large separation limit ||. The FCRLs they observe do 
not show spiral ordering of spins. 

We find that as v increases towards 1, the system 
changes (at non-zero g) smoothly from the TL phase 
characterised by (6,q)= (120°, 120°) to FCRLs with spi- 
ral orientational order. The FCRLs have different values 
of 9 and q, but their 9 + q value is always 240°, which 
is the value of 9 + q in the TL phase. As v is increased 
further there is another sharp transition from the FCRLs 
to the SL phase characterised by (9, q)=(90°, 180°). The 
parameter 9 + q shows a sharp jump from 240° to 270° 
at the FCRL-SL transition. Hence it is a convenient or- 
der parameter for this transition. The spin polarisation 
of the system varies smoothly with v and does not show 
an abrupt behaviour at any of the two transitions. The 
FCRL phases seem to smoothen out the jump in the spin 
polarisation, which was observed at the transition to the 
SL phase in previous work |@j. 

We consider skyrmionic excitations about the v = 1 
ground state. The number density of skyrmions is given 
by n s k y — ^-^n c , where n c is the carrier density. In 
our calculations we use the carrier density, n c — 1.5 x 
10 11 cmT 2 and change v by tilting the magnetic field. 
We localise the skyrmion centers on the lattice points 
of a FCRL shown in Fig. 1. 9 is the angle between the 
Bravais lattice vectors and e is the length of the Bravais 
lattice basis vectors, which vary with the filling factor. 
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We choose a skyrmion density (n s k y ) of 1 skyrmion per 
unit cell (i.e one skyrmion per lattice point). 

We measure all lengths in units of the magnetic length 

l c = \J~-^> and all energies in terms of the cyclotron en- 
ergy hto c , where lo c = (eB/m*c), m* is the effective mass 
of the electron. The area of the unit cell is fixed by the 
filling factor and is given by 
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The local spin polarisation which is represented by unit 
vector field n(x) is stereographically projected onto the 
complex plane by the transformation w — cot(#/2)e 1 ^, 
where 9 and <f> are the polar angles of the spin vector 
n(x). In the rest of the paper we will work with the 
planar spin variable w. 

The topological charge density is given by 



p(x) = — eijn.(9jn x djii) 
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The topological charge, Q(x) = J D d 2 x p(x) = 1 (J a 
denotes integration over the unit cell). 

The low energy, long wavelength excitations about the 
v = 1 ground state are accurately described by the 
NLSM. The NLSM energy functional has to be minimised 
for different filling factors to get the minimum energy 
configurations of this lattice. The NLSM energy func- 
tional with Zccman and Coulomb interactions is 



E = E, 
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The gradient or the spin exchange term proportional to 
/ d 2 x\ 9jn(x) 2 , (i = x, y), is calculated as 

T f 2 (d x wd x w + dywdyW) 



(1 + ww) 2 

where 7 = TTf^/f^- The gradient energy density is E grn 
/A , where A is the area of the unit cell. 

This term alone is the pure NLSM and it has scale 
invariant solutions ||. 

The Zeeman term is proportional to the z-componcnt 
of the total spin i.e. to ^ J n d 2 x ^ 1+ ^ (where n z is the 
z-component of n(x)). In our units the average number 
of electrons is ^ . The z-component of the total spin is 



(Total spin) z = — d 2 x— — 

2?r J a (1 + w 



ww) 



Therefore, 



E z = 9 *-\d* X 



(1 + ww) 



(5) 
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where g* = 3l i^ . The Zeeman energy density is there- 
fore E z /A, where A is the area of the unit cell. 

The Coulomb energy density term is a term of the form 
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where e* = (e 2 /Kl c )(l/huj c ) and A tot is the total area of 
the lattice. 

The Coulomb term arises because the electric charge 
density is proportional to the topological charge density. 
Since the topological charge density explicitly appears in 
the above expression the spin orientation gets automati- 
cally tied to the Coulomb energy The four dimensional 
integral in the Coulomb term can be converted to a sum 
over the reciprocal lattice. 
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where p(Gr) — J xea p(x)e lGn X , lies in the recipro- 
cal lattice and A is the area of the unit cell. 




Fig. 1 Face centered rectangular lattice with basis 
vectors: el=(ecos 8/2, — e sin 0/2),e2=(e cos #/2, e sin 6/2), 
nn,nnn,nnnn indicate nearest neighbour,next nearest neigh- 
bour and third nearest neighbour respectively. The box is a 
unit cell. 



We use the following ansatz to minimise the energy 
functional: 



{R} 



Ae 4qR 
R 



(9) 



where, A sets the scale for the skyrmion size, q.R is the 
spiral angle. The angle 8 characterises the lattice. This 
analytic ansatz exactly minimises the gradient term and 
gives scale invariant solutions Q . For analytic skyrmions 
the gradient term does not influence positional or spin- 
oricntational ordering. For overlapping skyrmions at low 
g, we still use the same ansatz since their size is large and 
the exponential damping term used in ref. [0 does not 
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have any significant effect. The size is determined by the 
competition between the Zeeman and Coulomb terms. 
We find that the Zeeman term prefers the SL phase and 
the Coulomb term prefers the TL phase. We expect to see 
FCRLs with spiral ordering when the transition occurs 
between the above phases. The NLSM energy functional 
has to be minimised with respect to the four variational 
parameters A, q(two components). 

To get an idea of the kind of positional and magnetic 
ordering expected in this system we look at a toy model 
, a spin system described by an XY model on the FCRL 
shown in the figure. The free energy of the system is : 
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We look at position and magnetic configurations which 
lie between the TL and SL phases, which are FCRLs. We 
consider only nearest neighbour (Ji), next nearest neigh- 
bour (J2) and third nearest neighbour (J3) couplings on 
the lattice. These configurations have 4 nearest neigh- 
bours and 6 next nearest neighbours (except at the square 
and triangle end points). We choose a spiral state with 
spiral angle q.R = q x R x +q y R y , where R x and R y are co- 
ordinates of a lattice point and q x and q y are the angles 
through which the spin turns for every jump to a lattice 
point. For the above configurations (except exactly at 
square and triangle end points), the free energy has the 
form : 



E = 2 Ji cos{q x ) 4- 2 Ji cos(q y 



2 J 2 cos(q x + q y ) 
-2J 3 cos(q x - q y ) (11) 



This free energy is minimised with respect to q x and q y . 
In the relevant parameter ranges, we find the solutions 
to be of the form (q x , q y ) — (q, q). This motivates the use 
of (q x , qy) = (q, q) in the exponential term of the ansatz 
in Eq. (§. 

With the above combination the number of variational 
parameters has reduced to three (#, A, q). The variational 
ansatz given above is used in the NLSM energy functional 
and all energies are calculated and minimised numerically 
(to 1 part in 10 5 ) with respect to these parameters . The 
filling factor v is varied for different values of the gyro- 
magnetic ratio g and the phase diagram for the system 
in the v — g/go plane is obtained (go is the physical value 
of the Lande g- factor of the system) . 

The accuracy of the routine was tested by integrating 
the topological charge density over a unit cell to get the 
topological charge, which should be 1 since there is one 
skyrmion per unit cell. An accuracy of 1 part in 10 6 was 
noted in the calculation of topological charge. 

Another test of accuracy can be done at g = 0. The 
only energy in the system is the Coulomb energy and this 
sets the length scale in the problem. It can then be shown 
that (9, q) should be independent of v and that A oc e. 
We find that the configuration which favors minimum 
Coulomb energy is the TL phase, (9, q)=(120°, 120°) and 
verify that A cx e. 
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Fig. 2 Phase diagram in the v — g/go plane 
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Fig. 3 Plots of 6 + q — v values at different g/go ■ The solid 
lines are the 9 + q values at the TL and SL, the dots indicate 
the numerical values obtained. Note the jump in the vicinity 
of the square transition. 



The minimum energy configurations were found for dif- 
ferent filling factors at various g values. For g = there 
is no transition from the TL phase as the filling factor is 
brought towards one. The TL-FCRLs transition starts 
appearing at non-zero g. The second transition to the 
SL phase, (6,q) = (90°, 180°), occurs very close to v = 1 
when g — 0.lg a . Both the transitions shifted away from 
v = 1 as the g value was increased and the intermediate 
region of FCRLs broadens. The phase diagram of the 
skyrmion lattice in the v — g/go plane is shown in Fig. 2. 
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The FCRLs appear as intermediate phases in regions of 
the phase diagram where the NLSM is a good description 
of the quantum Hall system (i.e low g and high v regions). 
The 9 + q — v diagram shown in Fig. 3. The 9 + q value 
is at 240° at the TL and FCRLs and is at 270° in the 
SL. The sharp jump seen in 9 + q occurs at all non-zero 
g values. Hence, we propose 9 + q as an order parameter 
characterising the FCRL-SL transition. 

The spin polarisation (sp) of the system of skyrmions 
at a particular v is 
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The FCRLs seem to smoothen out the jump in the spin 
polarisation that was observed in ref. ||. The value of 
total spin at g = O.lg was found to be 27 at v = 0.99. 
This means that many spins flip and explains the nega- 
tive value of spin polarisation . Direct comparison of this 
value with those obtained in low g experiments [ fb||l4| 
is not possible since we measure total spin at non-zero 
skyrmion density whereas the experiments measure spin 
for a thermally activated skyrmion-anti-skyrmion pair at 
v = 1. 




Fig. 4 Spin Polarisation at different g/go indicated in box 
brackets. 



Novel FCRL phases with spiral ordering such that 
9 + q = 240° appear in the low g limit (overlapping 
skyrmions). The TL-FCRL transition is smooth charac- 
terised by 9 + q = 240°. The second FCRL-SL transition 
is sharp and 9 + q jumps from 240° to 270°, making this 
a convenient order parameter for this transition. There 
are no spiral phases with 9 + q between these angles. 

Our calculations do not explain why FCRLs with these 
values of (9, q) occur but think that this is due to the 
hedgehog nature of the skyrmion. We have minimised 



the two-body energy functional of Timm et. al. and find 
that FCRLs with a general spiral ordering of spins do 
not occur. We confirm their phase diagram for a system 
of well separated skyrmions , which shows TL-FCRLs 
(Neel ordered)-SL transitions j^j. In their calculation, 
the Coulomb energy which favours the TL and the ex- 
change energy which favours the SL, compete to drive 
these transitions. We conclude that the spiral ordering 
we observe is not a two-body effect and occurs because 
of three-body or higher order effects. 

We thank Madan Rao and Surajit Sengupta for helpful 
discussions. 
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